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Abstract — By exploiting the freedom in defining the dual of spinors, we report an unexpected 
theoretical discovery of a quantum field theory of spin-half bosons. It fulfils Dirac’s 1969-70 ob- 
servation that “there must be boson variables connected with electrons”. The theory is local, 
Lorentz invariant, and has a positive-definite Hamiltonian. We formulate the unitarity-preserving 
scattering theory to accommodate the new dual and the associated adjoint. A model of Yukawa 
interaction with spin-half bosons and fermions of equal masses is studied to explicitly show uni- 


tarity. 
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There must be such boson variables connected with 
electrons. 


P. A. M. Dirac [1] 


Introduction. — With ninety-five percent of the exis- 
tence dark, in one way or another, its particle landscape 
cannot be a mere copy of one sort or the other of the 
Standard Model (SM) of high energy physics. This may 
appear as an audacious thought to many physicists. But 
the unexpected theoretical developments since 2004 [2-4], 
coupled with a lack of empirical signature for supersymme- 
try, suggests to consider the unexpected as a first-principle 
scenario for their darkness [4]. 

Recently, a new class of fundamental particles have been 
reported in [5,6]. While it is not widely known, Dirac him- 
self argued for electrons to have a bosonic partner [1]. 
Here, by exploiting a freedom in defining the dual of 
spinors (and the associated adjoint), we construct a the- 
ory of spin-half bosonic particles. These bosons cannot 
enter the SM doublets because of the statistics mismatch. 
The natural place for them is the dark-matter sector. As a 
by-product, they cancel the fermionic contributions to the 
cosmological constant if for every Dirac fermion of the SM 
there exists a here-reported bosonic partner of the same 
mass. 


(@)E-mail: dharam.v.ahluwalia@gmail.com (corresponding 
author) 


Quantum field theory of spin-half bosons. — In 
the SM, the Lagrangian density, at the level of kinemat- 
ics, is assumed rather than derived. The dynamics en- 
ters through the principle of local gauge invariance —the 
choice of symmetry group is then made on the basis of the 
hints hiding in the phenomenology. 


One exception to this general wisdom is to derive, rather 
than assume, both the quantum field and the Lagrangian 
density. This approach was initiated in 1964 by Steven 
Weinberg, and reached the textbook level in his mono- 
graph —The Quantum Theory of Fields [7]. 


His approach put the SM on a firm theoretical footing 
by unifying quantum mechanics with Poincaré space-time 
symmetries. An important departure from this approach 
occurred when mass dimension one fermions were first 
constructed [2,8]. Here we use the results from both ap- 
proaches to construct a bosonic quantum field of spin-half, 
and the associated Lagrangian density. 


Without assuming a Lagrangian density, we begin with 
the rest spinors for spin-half derived by Weinberg in the 
R ® L\j;=1/2 representation space [7] 


1 
&(0) = vm | 7 |, 
0 
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and 


0 

1 0 
£3(0) = vm o | $400) = vm 1 
—1 0 
These coincide with the eigenspinors of Q5 = yo, the fifth 
square root of the identity matrix in the ROL|;—=1/2 repre- 
sentation space [6]. The €;(p) for an arbitrary momentum 
are obtained by acting the R © L|;=1/2 boost operator, 


_ /E+m[l+g2 0 
D(p) = Im | O J= ae ? 


on €;(0): 


(3) 


&i(p) = D(p)é;(0), 


The crucial point of departure enters by defining a new 
dual [6,9] 


i=1,---,4. (4) 


&(p) = [P &i(p)] "yo (5) 
= telp)’ 1 1,2, (6) 
= [&:(p)] Yo: C= 3,4, 
where the covariant parity operator is defined as [3,10] 
pe m typ". (7) 


With the dual thus defined, the orthonormality relations 
and spin sums are found to be 


E(P) &j (p) = 2mdj,, (8) 
and 
XO &) Ep) = yp" + ml, (9) 
XO &(p)&(p) = —(yup" — ml) (10) 
1=3,4 


We now use the €;(p) as expansion coefficients for a 
spin-half quantum field, 


dp 1 


ao f OOJ | Do nPE)? 


i=1,2 


A U weiner (11) 


i=3,4 


This field is the same as the usual spin-half field of Dirac. 
The associated statistics follows only after its adjoint is 
defined. Instead of the Dirac adjoint, we introduce a new 
adjoint that exploits the new duals (6) 


3 
ata) © | Gor aE | Z, 10o 
pe holoe], (12) 
i=3,4 


where 


la:(p), a! (p+ = 5 (p — p') big, 


i 


lai(p), a; (p')]+ = [a; (p), a} (p')]+ = 0. (13) 


The upper and lower signs represent anti-commutator and 
commutator, respectively. We assume the same for b;(p) 
and bl (p). 

To establish the statistics we follow the procedure out- 
lined in [11]. Consider two space-like separated events x 


and y, and evaluate the amplitude for the particle, de- 


scribed by the a(x)-a(«) system, to propagate from x to y 
(with yo > zo). Since we are considering space-like sepa- 
rations, there exist a class of observers for whom 29 > yo. 
We call the former observers to be O and the latter to be 
O’. For the observers in O’, we calculate the amplitude 
for an antiparticle to propagate from y to x. Causality 
requires that these two amplitudes differ at most by a 
phase [11] 


Amp(xz > y, particle)|o = 


e? Amp(y — x, antiparticle)|o, (14) 


where 6 € R (see footnote '). An explicit calculation of 
the amplitudes in (14) using spin sums (9) and (10) while 
keeping track of the space-like nature of the (a — y) sep- 
aration, yields e? = 1. Thus, we must select the bosonic 
statistics. This is an unavoidable consequence of introduc- 
ing the new dual and adjoint in the a(x)-a(x) system. It 
evades the spin-statistics theorem [12,13]. The crevice in 
the theorem is the following: the new adjoint renders the 
interacting Hamiltonian non-Hermitian in a manner that 
still keeps the spectrum real. 

With the a(x)-a(x) system now fully specified, the 
Feynman-Dyson propagator for the new bosons, up to a 
factor 7, is simply the vacuum expectation value of the 
time order product 


=I 


( |¥la(x’)a(x)]| ) = ( |a(x’)a(w)| OE — t) 
+ (|a(w)a(x’)| ad — t’), 


with | ) as the vacuum state and Ẹ as the time-ordering 
operator. The two vacuum expectation values that appear 
in the right-hand side of the above expression evaluate, 
respectively, to 


(enata) = f (ga) ewe 


X (Yu p” + ml) (15) 
and 
7 1 
TOES f (=) eip (e—a) 
x [~ (7u p" — ml)], (16) 


1The “amplitudes” in (14) are dimensionful. It is the result of 
the fact that creation and destruction operators in the quantum field 
theory —for historical reasons— carry dimensions. 


24001-p2 


Spin-half bosons with mass dimension three-half: Evading the spin-statistics theorem 


where we have used the spin sums (9) and (10). Substitut- 
ing the two vacuum expectation values into the expression 
for ( |E[a(x’)a(x)]| ) and using the integral representation 
of the step function, 


1 co est 
t) = — 17 
O= eE (17) 
we obtain 
aN i —iq (a! —a 
(IENEN) = Geng f aeee 
Yuq” + ml 
— m. 1 
q-q— m? +ie (ta 


Consequently, the mass dimension of the new field —as 
referred to the a(x)-a(x) system— is three-half, with the 
Lagrangian density 

L(x) = a(x) (i770, — ml) a(x), (19) 


from which we obtain the canonical equal-time commuta- 
tors (thus establishing locality) 


[a(t, x), a(t, y)] =0= [p(t, £), p(t, y)] , (20) 
[a(t, æ), p(t, y)] = ið’ (æ — y)l, (21) 
where p(x) = ia(£)yo is the conjugate momentum. The 


free Hamiltonian evaluates to 


m= f oF 


x | J al(p)ailp)+ J bi(p)bl (p) 


i=1,2 i=3,4 


The free Hamiltonian is positive-definite if and only if the 
annihilation and creation operators satisfy bosonic statis- 
tics. 

The above expression can be rewritten as 


Ho = 25%(0) | ap |p|? + m? + [ëp |p|? + m? 
x | JO al (p)ai(p) + X of (p)bi(p) 
i=1,2 i=3,4 


Since i 
3 PEN 3 a 
6 (p) ~~ (27)? fe x exp(ip x), 


63(0) may be replaced by [1/(27)°] fda. With this ob- 
servation, the zero-point energy takes the form 


1 1 
Hac =4x ae) a 3p—V/|pi2+m?2. (22 
vac < aaa | 4 z [aps Ip? +m (22) 


In natural units A = 1 implies 27 = h, so Hyac represents 
a zero-point energy assignment of +4 /|p|? + m? to each 


unit-size cell (1/h*)d°x d3p in the sense of statistical me- 
chanics. The factor of 4 is consistent with the four degrees 
of freedom (when spinorial degrees of freedom associated 
with particles and antiparticles are accounted for). The 
second term in H shows that each of the four degrees of 
freedom contributes exactly the same energy to the field 
for a given momentum p. If for each of the SM Dirac 
fermions there exists a new boson of equal mass then the 
zero-point energies exactly cancel. This is a natural con- 
sequence of our formalism. 


Unitary evolution and scattering theory. — The 
bilinear invariants constructed from a(x) and a(x) that en- 
ter the interacting H are not Hermitian; see (26) and (27) 
below. Despite this, the interacting H still carries a real 
spectrum provided the norm of states is appropriately 
modified. 

To develop the formalism, we introduce 

Ht yH =H, (23) 
where ņ is a Hermitian operator, to be determined be- 
low. The work of Mostafazadeh then assures that H has a 
real spectrum provided that its eigenstates |y) have non- 
vanishing ņ-norms (Y|n|Y) # 0. To see this, we note 
that (23) yields 


0 = (b|(H'n — nF) iY) 
= (E* — E) (Vinly). 


A 0, 


(24) 


Therefore, when (yini) 
real [14-17]. 

To preserve unitarity, we construct 7 to have the follow- 
ing properties [16]: 


the eigenvalues are 


-1 


na;(p)n' = ailp), : 


nbi(p)n ~ = —b;(p). (25) 


That is, 7 commutes with a;(p), and anti-commutes with 
b;(p). Combined with the definition of the dual (5), this 
gives a relation between a(x) and a(z) 


na(x)n- = a(a). (26) 


Similarly, the bilinear invariant a(x)a(z) is related to its 
Hermitian conjugate by 


nla(x)a(a)|n~' = [a(x)a(a)]" (27) 
Since [ņ, Ho] = O, the free Hamiltonian is Hermitian. 


But in general, the full interacting Hamiltonian is non- 
Hermitian and does not commute with 7, so we re- 
quire (26) to ensure the reality of the spectra. The sim- 
plest contributions from spin-half bosons to H are of the 
form a(x)a(x) and a(x)7*a(z). 

An explicit solution for 7 that satisfies the above re- 
quirements may be constructed following [18]. To do that, 
we introduce ee 

n = eX, 


ae R, (28) 
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where 


x= [ 5 ol (p) 


1=3,4 


(29) 


Thus, 7 is a state-dependent relative phase between the 
particle-antiparticle sector. A short calculation reveals 
that 

— pia b! (p, o) 


nb! (p, a)n! (30) 


and helps us fix a. For eq. (25) to be satisfied, we must 
have a = m (modulo + integral multiples of 27). This 
gives us 7 = e’™X, which is Hermitian and unitary. To see 
that 1 is indeed Hermitian, notice nt = e~*™*', but for 
bosons y' = x; therefore ni = e~*™X. Now y acting on 
physical states |w) gives n where n is either zero, or an 
integer: xl) = n|w). Therefore, n'l) = e~'" |), while 
nlp) = et”). Thus while acting on physical states nt = 
n, since e~'™ = e’™. It is readily seen that 7?|wW) = |Y), so 
n? = 1. Consequently, while acting on |Y): nt = 77! = n. 

We now develop the scattering theory for spin-half 
bosons. Since for the states with time evolution governed 
by the generalized Hermitian Hamiltonian —of the class 
defined above in eq. (23)— their Hermitian inner product 
is not invariant under time translation 


(Ob liso = le ted) A hao BD 
It is necessary to introduce the 7-norm [14,16] 
(Al)n = (ln), (32) 
where 7 = eX. At t > 0, we find 
< (dln) liso = (OME — Hth) (83) 


Therefore, the 7-norm is invariant if and only if H satis- 
fies (23). 
Using the 7-norm, the S-matrix is given by 


SBa = (Bout |Qin)n 


= lim | tim (bol (r)n(7o)|a0), (34) 
where Q(T) = e*#7e~*#07 and we have taken the free 
Hamiltonian to be Hermitian. Therefore, 

SBa = lim „im olu (7, To) |Q0), (35) 
where 
U(r, To) = = nete —iH(T—To) eo tot | (36) 
In obtaining (36), we have used the identity 
cit Tye iHr =n, (37) 


which is equivalent to (23). 


Taking H = Ho + V and differentiating U (rT, To) with 
respect to T, we obtain 


[nV (T )n7 "| U(r, To), 


V(r) = et HoT VY e—*Hor | 


(38) 
(39) 


_d 
is U(r, To) = 


Using the initial condition U (To, To) = 7, the solution to 
U is given by 
z 

dr’ [nV 


U(r.) =9- 3 f (U, T) 


=n i — if dr'V(r') + 
TO 


A| w 


Taking V(t) = f d°xH(z), the Dyson series is 
sa =n +E A fatma dirn 
x (oha (x1): -- H(an)] |a0), (41) 


where nga = (boln|ao). 

Next, we address the issue of unitarity. For Hermitian 
Hamiltonians Ht = H, unitarity implies f dB! 5 Spa = 
ô(y—a) where f d@ integrates over momentum, sums over 
particle species and internal degrees of freedom. As for 
generalized Hermitian Hamiltonians, we have H#* = H so 
it suggests the following generalization to unitarity [19]: 


J d 88% SBa = 5(y — a), (42) 

where the matrix element of S¥, 76 is given by 
Sa = lim lim (yolU*(r,70)|G0), (43) 
U* (T, To) = 9 Ul (T, To)n. (44) 


To prove (42), we start with orthonormality and complete- 
ness relations associated with the free (and Hermitian) 
Hamiltonian Hp, and then implement evolution due to the 
full Hamiltonian H satisfying (23) as follows: 


J Aalao)(aol =, (Bojao) = 68-0). (45) 
Therefore, eq. (42) is satisfied when 
U*(r,79)U (7,7) = I, (46) 
for all r and 7. Rewriting U (r, 79) as 
U(r, To) = etHor p—tH" (7-10) p—iHoToy (47) 
we find 
U* (T, To)U (T, To) = 
eï HoTo [ee 70) ne iH" (7-70) | o iHoToy, (48) 
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Taking the inverse of (37), we get 


eH 70) le iHİ (r 7) = 9 1 (49) 
Since 7 = n7}, we obtain (46) as required. 
From (38) and (39), we obtain 
d 
-ir (1,79) = U# (T, To) VÝ (T). (50) 


Following the same procedure in deriving the Dyson series 
for Sga, we obtain 


Sfo met 5 J doatza 


x (yolnd [H (£n): H'(21)] |Bo). 


Normalizing the S-matrix as 


Sga = Nga 7 2riMgað (pg = Pa); (51) 
Si, = nye + 2miM¥,5" (py — pp), (52) 
eq. (42) becomes 
i fas [nya Maad* (pe — Pa) 1; MŽngað (pa = ps)| = 
Qn J dp [5* (we — p,)54(pp — pa) MŽ; Mpa] (53) 


Setting y = a, the term 64(pg — pa) cancels from both 
sides of (53) so we obtain 


i fas (nap Mpa -M #onpa) = = 


Qn / dB [5*(we ~ pa) M# Mpa (54) 
Equations. (53) and (54) can be seen as a generalization 
to the optical theorem. In the subsequent section, we will 
demonstrate that (54) is satisfied for a model of Yukawa 
interaction involving spin-half bosons. 

To compute physical observables, we first need to define 
the transition probability for the process a > 8. Since the 
unitarity of the S-matrix is defined using the generalized 
Hermitian adjoint #, the transition probability should be 
appropriately modified. In a finite volume V with duration 
T, we define the transition probability to be 

P(a > p) =VT(2n)* (GapMSsMpa)- (55) 

Because M#,M, ga is not guaranteed to be positive- 
definite, we multiply it by a constant phase pga E€ R to 
ensure that P(a —> 8) > 0. From the definition of 7 
and the Dyson series of Sga and S7 we find that the 
Mz, = eit(natns) MI g where na and ng are the number 
of spin-half anti-bosons present in states |ao) and |80), 
respectively. Therefore, by choosing pag = e77 ate), 
eq. (55) reduces to the standard transition probability and 
the physical observables defined using Hermitian conjuga- 
tion can be applied without difficulty. 


Phenomenologies. We now consider two spin-half 
bosonic and fermionic fields of equal mass coupled to a 
real scalar field ¢(x) 


=o [ae w(x 


where y(x) is a fermionic field whose expansion coefficients 
are given by (1) and (2). By the construction of 7, it 
commutes with the fermionic and scalar fields 


(x) + a(x)a(x)| (x), (56) 


(57) 


Therefore, the Yukawa interaction (56) is a generalized 
Hermitian operator n~!V1(t)n = V(t). 

We now demonstrate that the Yukawa interaction sat- 
isfies (42), the generalized unitarity relation. Taking the 
initial state to be a single scalar boson, the matrix element 
of 7 is Nea = nge = 0(6 — $) so (54) becomes 


i(Mse — MŽ) )= 2r fas [ori (pe — po) M$,Mpo| - (58) 


We will evaluate both sides of (58) independently. Starting 
with the left-hand side, the matrix elements Mgq and Mz, 
are obtained by evaluating the amplitude for ¢ > ¢’ and 
subsequently setting ¢’ = ¢. Using the solution of 7, we 
obtain 7|¢) = |¢) where |¢) denotes a single scalar boson 
state. Setting ¢ = ¢’, we find Mi, = MŽ. Therefore, we 
may rewrite (58) as 


Im(Mgo) = =r | a8 [3" (po - pe)M$yMao]. (69) 
Here, Im(M gq) is proportional to the imaginary part of 
the loop corrections to the scalar propagator [20]. Due 
to the difference in statistics and the fact that the spin- 
half bosons and fermions have equal mass, the bosonic and 
fermionic loop corrections identically cancel to all orders 
in perturbation so Im(Mg¢¢) = 0. Next, we evaluate the 
right-hand side of (58). It expands to 


fæl [5*( (pp — Pa) M¥,Mp¢| = 
4 | 
i 2 Sru Mayet Pg + Pv — pe) 


OSSA aT ô (Pmt Pa — Po), 


spins 


(60) 


where dy = d’pyd’pz and da = d?padp— The first 
a 
and second term of (60) comes from the fermionic and 
bosonic interactions, respectively. Acting 7 on the relevant 
states in (60), we obtain 
nlo) =l); nldv) = Yy), njaa) = —|aa), 
where Yy) and | aa) denote the pair of spin-half fermion— 


anti-fermion and boson—anti-boson states, respectively. 
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The actions of 7 on |) and jaa) have opposite signs 
which induces a relative phase in (60). In effect, we find 


Data availability statement: No new data were created 
or analysed in this study. 
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